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TETHERED SUBSATELLITE STUDY 
I. INTRODUCTION 

Mathias P. L. Siebel 


Various concepts of using tethers in space have been studied previously; 
e. g. , an earty concept for Skylab had the solar observatory (Apollo Telescope 
Mount) tethered (rather than hard-docked) to the S-IVB stage, and other tether 
studies were performed in connection with astronaut rescue concepts. The 
present in-house study was performed to gain an understanding of a report 
received in September 1974 from the Smithsonian Astrophysical Observatory 
(SAO ) entitled Shuttle-Borne "Skyhook": A New' Tool for Low-Orbital-Altitude 
Research, by G. Colombo, E. M. Gaposchkin, M. D. Grossi, and G. C, 
Weiffenbach. 

In their report SAO suggests a tether of approximately 100 km length 
deployed from the orbiter and carrying a subsatellite to an altitude of perhaps 
100 km below the shuttle and 100 to 150 km above the Earth (Fig. I-l). The 
subsatellite would then be used as a platform on which various kinds of experi- 
mental apparatus would be carried for upper atmospheric measurements, high 
resolution gravity gradient measurements, and low frequency radio physics 
experiments in the lower magnetosphere where the tether wire itself is used as 
a dipole. 

An in-house study performed by C , C. Rupp (see Bibliography) verified 
Certain parameters and I’esutts of the SAO report. A further, more detailed 
in-house study was performed during the period from February 1975 to August 
1975. The following are the principal aspects and the responsible investigators 
to the tether-subsatellite system that were considered: 

Dynamic Analysis of a letliered Subsatellite Mario H. Rheinfurth 

and Zachary J. Galaboff 

Investigation of a Tethered Subsatellite Ralph R. Kissel 

Control Law 


Aerodynamics 


Kenneth D. Johnston 



Thermal Studies William P. Baker 

Tethered Subsatellite Communications J. A. Dimliin 

System Design 

Reports on these individual efforts form the body of this repor*:. 

It will be seen that while many detailed questions remain imresolved, no 
fundamental effect has been discovered that makes the subsatellite concept 
unfeasible. More detailed systems and optimization studies are recommended. 
Further, it should be mentioned that in the course of the study, interest in using 
the tether was expressed by a number of members of the scientific and applica- 
tions o^’iented commimities. In addition to the applications originally suggested 
by SAO, other uses for the tether have been found both in the Atmospheric 
Magnetospheric Plasma in Space (AMPS) area (generation of Alfven waves and 
study of plasma waives) and in more general areas. The more general applica- 
tions proposed include cargo transfer, space station stabilization, the use of 
tethers as aids in the erection of large area space structures, satellite retrieval, 
power generation, and others. 

It is recommended that further studies be performed to define common 
and specific application function features of tethered subsatellite systems. Such 
studies should culminate in an experiment on an early shuttle flight. 
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Figure 1-1. Tethered sutasatellite. 





1 1. DYNAMICS ANALYSIS OF A TETHERED SUBSATELLITE 

Mario H. Rheinfurth and Zachary J. Galaboff 


A. Summary 

A dynamic analysis was performed to study the feasibility of the tethered 
subsatellite concept with emphasis on its dynamic response and stability charac- 
teristics. The differential equations were derived that describe the three- 
dimensional motion of the subsatellite and the flexible tether that connects it 
with the main satellite. Because of limitations in manpower and resources, 
only a simplified model could be simulated on the computer. However, a serious 
effort was made to include the most significant factors in the analysis, factors 
that were deemed necessary to prove the dynamic feasibility of the concept. 
Analytical investigations and computer simulations performed thus far did not 
expose dynamic characteristics that would preclude a satisfactory deployment 
and stationkeeping of the tethered subsatellite. However, a serious dynamic 
problem was encountered during the attempt to retrieve the subsatellite. This 
problem occurs in the vicinity of the subsatellite and is characterized by a 
strong tendency of the subsatellite to sling aroimd the main satellite. This 
phenomenon and the lateral tether flexibility, which was not included in the 
present analysis, are still critical issues and require further study. 


B. Introduction 

The objective of this study is to conduct a feasibility analysis of a sub- 
satellite that is attached to an orbiting spacecraft by a tether. For this puiTose, 
a computer program was developed that allows the prediction of the dynamic 
behavior of the tethered subsatellite during its deployrr.ent, stationkeeping, and 
retrieval. The dynamic analysis is kept sufficiently general to accommodate a 
rather wide variety of system parameters. Because the equations of motion 
are fairly complex for the general case of three-dimensional motion and tether 
flexibility, they are written in concise vector-dyadic and matrix form using 
Lagrangian mechanics. The derivation is carried to a point from which the 
interested reader can readily proceed towards a detailed scalar formulation. 

The computer simulation itself was performed under certain restrictive assiunp- 
tions dictated by the nature of the study and by the limitations of manpower and 
resources. A detailed formulation of the corresponding simplified equations of 
motion is provided. In selecting the simplifying assumptions, particular 



attention was given to the identification and deletion of dynamic effects tiiat can 
be safely neglected. A point-by-polnt discussion of these effects is included. 
Consideration is also given to effects that should be included in a more refined 
analysis. The necessary extension and augmentations of the equations of motion 
can be obtained without difficulty from the vector-dyadic formulation presented 
in this section. 


C. Lagrange's Equations 

The dynamic analysis of complex systems is greatly facilitated by an 
approach attributed to Lagrange and generally referred to as anal 5 rtical 
mechanics. This approach is general and systematic in nature and is readily 
adaptable to modifications and refinements of the mathematical model that is 
used to describe the dynamical system. It is obviously beyond the scope of this 
report to derive the equations of motion in an expository manner. A ceiTain 
familiarity of the reader with the concepts of analytical mechanics will, there- 
fore, be required for the understanding of this derivation. However, this 
familiarity is not necessary for the application of the equations of motion as 
they appear in their final form. On the other hand, tl)e treatment is detailed 
enough that the cognizant reader should be able to supply potentially desirable 
addenda to the existing mathematical model without great difficulty. 

In their classical form, Lagi'ange's equations are based on an inertial 
reference frame and employ generalized coordinates. Accordingly, they appear 
in the form of a matrix equation 


_d /8T\ 

dt J 


8T 

hq 


+ X = Q 


(Il-l) 


where T = T(£,£) is the kinetic energy, expressed as a function of tlie general- 
ized coordinate vector q and its time derivative q. The term X represents 
the Lagrange multiplier vector, and C represents tiie constraint matrix 
appearing in tlie matrix equation for the constraints liaving the Pfaffian form; 


‘I> = Cq - b - 0 


( 11 - 2 ) 


II-2 



The generalized force vector Q on the right side is determined via the prin- 
ciple of Niirtual work; 


6V/ = ^ 5 ^ 


(II-3) 


This states that the viiti’ai work 6W of the applied generalized forces ^ is 
zero for virtual displacements 5^ which are consistent with the constraints 
imposed upon the system. 

Equation (II-l) represents a set of second-order ordinary differential 
equations. The solution of these is necessary and sufficient to establish the 
complete dynamical behavior of the system as a function of time. However, it 
results that for more complex dynamic configurations, this form of the equations 
of motion is too complicated. A less complicated form of these equations is 
obtained by introducing one or more noninertial reference frames and expressing 
the kinetic energy in terms of nonholonomic velocities. These velocities are 
also referred to as derivatives of quasi-coordinates. The corresponding trans- 
formation from generalized velocities to nonholonomic velocities is given by 
Euler's kinematical equations 


0 = A(q) a 


(II-4) 


with A(q) being the appropriate transformation matrix. The nonholonomic 
velocity vector H is, in general, composed of both linear and angular velocity 
components. Introducing equation (11-4) into equation (II-l) yields Lagrange's 
equation in a quasi-coo nlinate form as 


_d 

dt 



+ 


((A - J) A"‘ 


T 9T 

0g 


(A-') 


T ^ 

9a 


T 

B \ 




where B = CA"‘ 


and where J 


is tlie Jacobian matrix 


(II-5) 


J = ao/0£ 


( 11 - 6 ) 


n-3 



Judged by their outward appearance, the transformed equations (11-5) seem to 
be more complicated than tlieir classical coimterpart given in equation (ll-l). 
However, th ir intrinsic simplicity will be obvious when the detailed steps of 
introducing the above-mentioned noninertial I'efei'ence frames are carried out. 
Since these steps are rather lengthy — but straightforward — they will not be 
repeated at this point. They will produce the result that for each Introduced 
noninertial reference frame, tlie equations of motion can be partitioned into 
three distinct sets. These sets can be physically interpreted as rigid-body 
translation, rigid-body rotation, and subsystem flexibility. Using vector-dyadic 
notation, this result can be stated in the following fonn: 

Rigid-Body Translation 


4 + wx^+r *X=F 

dtlSv/ — 8v V - — 


Rigid-Body Rotation 


_d 

dt 


/8T\ T 

I — I + wx — 

y0w y ~ 

Subsystem Flexibility 


•7^ + r • X = L 

w 0) ~ — 


(11-7) 


( 11 - 8 ) 


~ l~\ - ~ + T • 
dt \9q / q 

where F , F , and F are the pertinent constraint dyadics as determined 
V w q 

by the properly transfonned constraint equation (11-2) . 

The elimination of the unlmown Lagrange multipliers is usually performed 
by a separate computational routine. A frequently used one is outlined in the 
appendix. 


A = Q 


(11-9) 


n-4 
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D. Coordinate Systems 

1. ORBITAL REFERENCE FRAME 

The motion of the tethered subsatellite will be described relative to an 
orbital reference frame that is fixed in the main body. Since it is known in 
advance that the short-term dynamic effects of a tethered subsatellite on its 
main satellite are by engineering design extremely small, the dynamic analysis 
can be greatly simplified by assuming that the reference frame is traveling 
along a predetermined orbit. This orbit can be separately calculated by a 
trajectory analysis. Except for rather low orbits, this analysis requires only 
the consideration of gravitational forces. For a spherical gravitational field, 
the orbits uill be Kepleiian. The orbital reference frame is aligned such that 
the x^-axis points in the direction of the inertial velocity vector, the z^-axis 
toward the center of the Earth, and the y^-axis completes the right-hand triad. 

It is a noninertial coordinate system. In the particular case of a circular orbit, 
the x^-axis will be parallel to the inertial velocity vector. The long-term . 
dynamic effects of the tethered subsatellite on its main satellite can be calculated 
with sufficient accuracy in a separate analysis after tire completion of the 
dynamic analysis described herein. 

2. QUASI-COORDINATE REFERENCE FRAME 

To derive the equations of motion using the quasi-coordinate formulation 
of Lagrange's equations, a noninertial quasi-coordinate reference frame is 
introduced. The axes of this second reference frame will be labeled x, y, z. 

Its origin is chosen to coincide with the orbital reference frame. It is 
advantageous to align the axes of the quasi-coordinate frame parallel to the 
principal axes of the undefoimed tether/subsatellite body. The orientation of 
the quasi-coordinate reference frame relative to the orbital reference frame 
can be defined by three Euler angles. Initially, the two frames are coincident. 

A series of three rotations in the proper sequence defines the relative orienta- 
tion. The three rotations are: 

a, A positive rotation iji about the z-axis. 

b, A positive rotation 0 about the y- axis, 

c, A positive rotation ^ about the x-axis. 

Frequently, these Euler angles are given the following names; the yaw (heading) 
angle ip , the pitch (attitude) angle 0, and the roll (banl^) angle (p. Figure II-l 



illustrates the relationship between the two reference frames. For the salce of 
visibility, only the last two angles are shown in the figure. The angrilar velo- 
city w of the quasi-coordinate frame relative to the orbital reference frame 
can be e;^^essed in terms of the time derivatives of these Euler angles through 
the well-lnioWn Euler kinematical equations. This relationship is given as; 


= 0 - ip sino 


CO = 0 cos 0 + d cos 0 sin ' 

y 


o> = 0 cos 0 cos <t> ~ 0 sin < 
z 


( 11 - 10 ) 


The Eulerian angles define the orientation of a body by the minimum number of 
independent coordinates. Their primary attractiveness stems from the fact 
that tliey permit a simple geometrical interpretation of the attitude motion of 
a body. However, they suffer from an intrinsic singularity condition. Another 
computational drawback of the Euler angle method is the nonlinearity of the 
differential equation (II-IO) which has to be solved to obtain the Euler angles. 
Both of these computational disadvmtages can be avoided by defining the orienta- 
tion of a body in terms of direction cosines or quaternions. Since one of die 
primary purposes of this study was to gain physical insight into the dynamic 
characteristics of the tethered subsatellite, the equations of motion were, 
nevertheless, derived in terms of Euler angles. 


E. Kinetic Energy 


The kinetic enei-gy of the dynamical system relative to inertial space is 
defined as 


T 


1 

2 


/tZo 


+ 


(Qx n) 


• 2 

+ RJ 6m 


( 11 - 11 ) 


n-G 



where 


Vq = inertial velocity of the origin of the quasi-coordinate frame 

n = angular velocity of tlie quasi-coordinate frame relative to inertial 
space 

^ = position vector from the origin to a mass particle 
E = velocity of a mass particle relative to the quasi-coordinate frame. 


The integration extends over the total tethei’ed subsatellite system. 

Because the origin of the orbital reference frame is assumed to be in a 
state of free-fall, its velocity can be set equal to zero ( Vq = 0) . A formal 
proof of tliis statement will not be presented. The angular velocity n of the 
quasi-coordinate frame can be expressed as 


^ + w , (11-12) 

where 

X = angular velocity of orbital reference frame relative to inertial 
space 

w = angular velocity of quasi-cooi’dinate frame relative to orbital 
reference frame. 

The position vector R can be expressed in thre§ components as 

R = f 0 + I + w (11-13) 
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where 


r = position vector of a mass particle originating from ^ 

w = elastic displacement vector of a mass particle frdih undeformed 
state. 


The generic vector is introduced as an extraneous coordinate for the purpose 
of determining internal x’eaetion forces at points of int^lfest. The elastic dis- 
placement vector will be described in terms of the n i%i^ modes (eigenfurt.ctions) 
of the dynamical system; i. e. , r ^ ' 'T; , , : ' 




where represents a three-dimensional normal mode and I’epre- 

sents its associated generalized cooi'dinate. The normal modes are calculated 
by a separate structural analysis. Using normal mode information greatly 
facilitates the formulation and solution of the equations of motion. Although it 
is realized that such a modal analysis is theoretically only applicable to linear 
time-invariant structures, past experience has demonstrated that it can also be 
used for systems whose parameters change only slowly with time. The tethered 
subsatellite system falls into this category because tire deployment and retrieval 
velocities have to be kept sufficiently small to prevent undesirable nonlinear 
response phenomena and dynamic instability of the tether because of excessive 
Coriolis effects. 


F. Equations of Motion 


Since the objective of this study was to explore the feasibility of the 
tethered subsatellite concept, the analysis was based on a rather simplified 
dynamic model of the system. However, some of the simplifications had to be 
introduced because of limitations in time and resources. The most significant 
simplifications in this latter category were the omission of laterjil tether 
dynamics and the use of only one normal mode for tlie elastic expansion of the 
tether. Torsional rigidity of tlie tether was also neglected. However, its 
effect will be voi^ small. 
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In general, the elastic defoimation vector of the tether cfln be expi'essed 
in modal form as follows: 

m 

w(r,t) = [S X.(x,y,z) 5 .(t)i et + [2 Y.(x,y,z) I £2 


+ [2 Zj^(x,y,z) t,|^(t)Je3 


( 11 - 16 ) 


where e^, 62, and 63 are unit vectors along the axes of the quasi-coordinate 
reference frame. The terms Xj(x,y,z) and Y.(x,y,z) represent the lateral 

normal tether modes v;ith ^^(t) and hj(t) beii.^^ their corresponding generalized 

coordinates. The term Zj^(x,y,z) represents the longitudinal normal tether 

modes and t,j^(t) ths^^orresponding genei’alized coordinates. With the above- 

mb^iloned simplifications, the deformation vector reduces to 


w(r,t) = Zi(z)Ct(t)e3 


( 11 - 16 ) 


The first longitudinal mode was appx’oximated by 


Zi(z) 



(n-lV) /g 


This approximation was considered sufficiently accurate for the feasibility study 
of tlie tetliered subsatellite concept. 

The tether itself was assumed to have a constant diameter and imiform 
mass distribution per unit length. For simplicity, the subsatellite was given a 
spherical shape and mass distribution. Because tether torsion was not taken 
into consideration, the rotational motion of tlie subsatellite abovit the tether axis 
was set equal to zero, Thei’efore, the motion of the subsatellite can be described 
In tenns of omy two of tlie three Euler angles. The two Euler angles chosen 
were the pitch angle 0 , which describes the in-plane motion of the subsatellite, 
and the roll angle <(i , which describes its out-of-plane motion. The equations 
of motion can tlien be expi’essed in tlie following form: 


II -9 


In-Plane Motion 


I^M + ^ cos^0 + ij 0 + 2 ^ ^ dj L( (5 - y 

2^M + L^ (y - o) 0 sinc^ cos-0 = 


y) cos'0 


(11-18) 


Out-of-Plane Motion 


|Mtf)l,>n]y . 2[(wf)«' • ji 


L 0 


(m.s) 


(y - O)^ cos 0 sin0 = Q 


<P 


Stretch Equation 


(-f)f 


I M + ^ f 


(11-19) 


I.jj/*" ' (y - fl)^ cosVj = Q, 


( 11 - 20 ) 


The ation for tlie tension in tlic tether is obtained by differentiatin"' the 
Lagrange equation (II -9) with respect to the extraneous variable If), For the 
design of the tether control system, it is nocessai-y to know the tether tension 
at the deployment reel/tether intei’facc. Tliis is obtained by setting the 
extraneous valuable 


ijo - -^0 S) 


( 11 - 21 ) 


and the corresponding constraint condition to 


e. = 


( 11 - 22 ) 
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where represents a prespecified time function that is detennined by the 

particular tether control law. The Lagrange multiplier A appearing in equation 
(n-9) yields the tension in the tether directly. After all necessary differentia- 
tion steps are performed, the tension equation is obtained by setting io ~ 
ig = i , and io = f . This yields 

(M + m)i + + y^L[^^ + (y - O)^ cos^0l + T = 

( 11 - 23 ) 

where T is the tether tension. 


G. Generalized Forces 


The generalized forces arise from several different sources. The most 
significant of these will be discussed in this subsection. Others can be added 
later if needed. As previously mentioned, the generalized forces are calculated 
by the virtual work done by the external forces through virtual displacements. 
The virtual work can be conveniently expressed in terms of the coordinates 
(^> Yo’ H) Oie orbital reference frame and the corresponding components 
of the external forces. The desired generalized forces are then obtained by a 
subsequent transformation to the generalized coordinates. Thus, the vix’tual 
work is 


6W=F6 +F6 +F6 

^0 ^0 ^0 ^0 '"O =^0 


( 11 - 24 ) 


The position of a mass particle of the tether can be defined in terms of the 
generalized coordinates of the dynamical system as 


Xo = 

(f„ + z + $(z) rj 

sin 0 cos 

yo = 

- [to + z + <I>(z) ^ 

1 sin^i 

zo = 

[fo + z + $(z) C 1 

COS0 COS0 


( 11 - 25 ) 
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where $(z) is the approximation used for the first longitudinal mode of the 
tether. For the calculation of the generalized forces due to a virtual displace- 
ment of the subsatellite, equation (n-25) is used setting z = i and $(i) = 1. 

The virtual displacements can be expressed in terms of the generalized 
coordinates by using the following differential relationships: 


6 = [fg + z + $(z) cos 0 cos <l>dO - [fg + z + $(z) C 1 sin 0 sm f 6^ 

^0 

+ sin 0 cos g + $(z) sin 0 cos f 6^ 


6 = - [fg + z + 4>(z) cos <p 6<p - sin ^ 6^g - <h(z) sin ^ 6^ 

^0 

(11-26) 

6 = - [ f g + z + $(z) 4 1 sin 0 cos ^ 6o 

- [fg + z + $(z) ^ J COS 0 sin (p dtp 


+ cos 0 cos <p 6fg + 4>(z) COS 0 cos <p 6^ 


Inserting equation (II-2C) into equation (11-24) and collecting the appropriate 
terms yields the virtual woi’k as 


6W = do + 60 + 6^ + 6fp . (II-27) 


1. FIRST-ORDER GRAVITY GRADIENT FIELD 

Since the origin of the orbital reference frame moves along a free-fall 
trajectory, the only gravitational forces acting on tlic tethered subsatellite arise 
from the gravity gradient field. The gravity gradient force terms arc obtained 
by a Taylor-series expansion of the gravity field about the free-fall trajectory. 
The first-order terms of this series are well knovm. Applying these terms to 
a mass pai-ticle of size dm results in the following; 
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(1) 9 

dF ^ = - ojg 5<o dm 


dF = - coj yo dm 
^0 


dF = 2 Ct)j Zo dm 
^0 


where 


OJo 


gP^E 

R3 


(11-28) 


It should be noted that these expressions are not restricted to circular orbits, 
but are valid for general orbits followed by the orbital reference frame. Sum- 
ming up the forces over all mass particles of the dynamical system yields the 
first-order gravity gradient terms as: 


0 

g(Gl) 




wj L* sin 0 cos 9 cos^0 


Wp iJ cos^ 0 sin <p cos <l> 


L(3 cos^ 0 cos^ 0-1) 


Q 


(Gl) 




L(3 cos* 0 cos* 0-1) 


(11-29) 
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2. SECOND -ORDER GRA\aTY GRADIENT FIELD 

The accuracy of the dynamic model can be improved by including higher 
order terms of the gravity gradient field in the analysis. For most practical 
space structures, however, it will not be necessary to go beyond second-order 
terms. These are less known than the first-order terms, but their derivation 
is straightfoi'Nvard. Applying these second-order tenns to a mass particle of 
size dm yields the following forces: 


dF 


( 2 ) 


- 3 CO 


2 

0 


Rfl 


dm 



- 3 w 


2 

0 ■ 


yn ^fl 


dm 


dF 





(11-30) 


In comparison with the first-order gravity gradient terms , it is seen that they 
are very small. Even for distances up to 100 km from the origin of the orbital 
reference frame, their contribution is on the order of a few percent. They were, 
therefore, not considered in the analysis. 

The generalized forces associated with these forces can again be obtained 
by summing over the whole dynamical system and applying equations (11-25) and 
(11-26) . 

3. ROTATING ATMOSPHERE 

For tethered subsatellites moving in orbits of extremely low altitudes 
( 100 to 200 km) , the aerodynamic disturbances arising from the atmosphere of 
the Earth dominate those arising from the gravity gradient field and decisively 
influence the dynamic response of the system. To limit the size of the computer 
program, the variation of the atmospheric density with altitude was assumed to 
be exponential. A least-square curve fit for such an "exponential" atmosphere 
is shown in Figure 11-2. For the analysis, it was of some advantage to define 
the exponential density vax'iation in the form 
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(11-31) 


bli h/H 

p = Po e = Po e 

where 

Pq = reference density 
H = ^ = reference height 

h = distance from main satellite vertically down. 


The data are based on the 1962 U. S. Standard Atmosphere. Only a single scale 
height H was used in the analysis, which corresponded to an altitude range of 
80 < h < 100 km. This was considered the range of major concern. As a con- 
sequence, the atmospheric density in the vicinity of the main satellite was much 
smaller than the actual one. Therefore, the dynamic response of the tethered 
subsatellite during its initial deployment phase was not accurately modeled. 
However, this was not considered to be critical, because it was found that the 
initial dynamic transients damp out rather quicldy during the deployment phase. 
In fact, an increase in drag on the subsatellite will benefit the deployment initia- 
tion. For the tether, the aerodynamic forces were defined as; 

a. Normal Force 


D =-pCAVV 
-N 2 ^ N N N -N 


b. Axial Force 


D. = -XPC A V\' 
—A 2 AAA ~A 


(11-32) 


(11-33) 


where V and V . represent the normal and axial relative wind velocities. 

— N —A 

The normal and axial force coefficients and are usually given as 

functions of angle of attack and other aerodynamic parameters. In the present 
analysis, they were assumed to be constant. 
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For the spherical subsatellite, the aerodynamic force was defined as 


D 



A V V 


(11-34) 


with V being the relative wind velocity and the aerodynamic drag coeffi- 
cient. The latter was assumed to be constant. 

The relative wind velocity was determined with the assumption that the 
atmosphere rotates with the Earth. Resolving the relative wind vector in com- 
ponents along the axes of the orbital reference frame yields the following 
results: 


V = - Ro (wo - cos i) 
^0 


V = Rfl (S2o sin i cos (3) 
^0 


(11-35) 


It can be easily verified that the motion of the tetliered subsatellite relative to 
the orbital reference frame is negligible when compared to the total relative 
wind velocity. Therefore, the velocity terms in equations (II-32) , (II-33) , and 
(11-34) can be directly calculated from equation (11-35). The detailed steps of 
this calciUation will not be given. Likewise, the calculation of the generalized 
aerodynamic forces acting on the system via the principle of virtual work will 
not be presented. The final x’esult gives the generalized aerodynamic forces as 

q(A) ^ p ^bh, cos » 

(II-3G) 
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^ '""["d * c;-) rp. Y^")] 


with the following definitions: 


V = V cos 0 

X Xq 


V = V sin 0 sin 6 + V cos 

y y„ 


(11-36) 

( Concluded) 


V = V sin 0 cos <t> - V sin <b 
^ ^0 ^0 


and 


,(M) 




^0 ^0 


V, 


■ 

N V X y 


V = V 
A z 


Also, 


1 bL cos 0 cos ' 
^0 " 2 ® 
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bL cos 0 cos (t> 
e - 1 

b cos 0 cos <j> 


Ml = 

bL cos 0 cos 0 / , . „ . , 

_ e ( bL cos <1 cos 0 - 1 ) + 1 

(b cos 0 cos 0 ) ^ 
h 

hfl = y [1 - cos (y+ Yo)] 


The last of these definitions is to be used for elliptical orbits of the main 

satellite where h represents the altitude variation because of the eccentricity 
e 

of the orbit. The angle y represents the true anomaly of the main satellite, 
and the angle Yo provides a convenient definition for initiating various maneu- 
vers of the tethered subsatellite. 

The numerical values for the aerodynamic coefficients were assumed to 
be constant at the following levels: 

Subsatellite — = 1.0 

Tether =2.2 
N 



4. ORBITAL ECCENTRICITY 

The effect of orbital eccentricity on the dynamic response of a tethered 
subsatellite is twofold; one is a geometric effect, and the other is a dynamical 
effect. The geometric effect is due to the altitude variation of the tethered sub- 
satellite which is caused by the eccentric orbit of the main satellite. As a 
consequence, the tethei’/subsatellite system is subject to unsteady aerodynamic 
distui’bances. To incor{}orate this effect into the analysis, it was assumed that 
the main satellite moves in a Kepleiian oi-bit whose geometry is defined as 
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R (1 + e) 

—JL 1 

1 + e cos y 


(11-37) 


where is the perigee distant, e is the eccentricity of tlie orbit, and y is 
the true anomaly. 

Since the eccentricity e will be small (e « 1) for all orbits of practical 
interest, it is possible to expand equation (11-37) in a Taylor series as follows: 

R = (1 ~ e y+ e^ cos^ y - . ..) . (11-38) 


Taking only the first-order teirm in e , the altitude variation becomes tlien 


Ah = R - R = e R (1 - cos y) 
P . P ' 


(11-39) 


From this it follows that, for small eccentricities, the altitude variation of the 
tethered subsatellite is approximately sinusoidal. The relationship beRveen the 
maximum altitude change and the orbital eccentricity can be readily obtained 
from equation (11-39) as 


e 



(11-40) 


It is instructive to consider a numerical example. Assuming an altitude of 
200 km for the main satellite (Kp= 6^70 Ion) and a maximum altitude variation 

(h^= 20 km) of the subsatellite leads to an orbital eccentricity of e= 1. 5 x 10"®. 

This result confirms the validity of the above Taylor-series expansion. 

Having demonstrated the smallness of the eccentricity, it is now possible 
to assess the dynamic effect of the eccentricity using a series expansion express- 
ing the true anomaly y in terms of the eccentricity e and the mean anomaly 
M . This relationship is given by 
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y = M + ^2 e - I- sin M + ^ 4 2M + 


where 


^go K 


M = 


t = nt 


and 


n = 27 t/T 


(11-41) 


with a representing the semimajor axis and T the orbital period. Neglecting 
terms higher than first order in equation (11-41) yields the angular velocity of 
the main satellite in its orbit as: 


y = n (1 + 2 c cos nt) (11-42) 

From this it follows that the angular velocity of an eccentric orbit is sinusoidal 
for small values of its eccentricity. However, it is seen that the sinusoidal 
fluctuations are extremely small. For most practical orbits, these fluctuations 
are less than 1 percent of the mean angular velocity. It is, therefore, permiss- 
ible to replace the time-varying angular velocity y by its mean angular velocity 
n in the equations of motion. 

The angular acceleration of the eccentric orbit is obtained by differen- 
tiating equation (11-42) with respect to time. This yields 


y = -2 e n^ sin nt 


(11-43) 


The inertial forces arising from this angular acceleration are, therefore, piu- 
portional to c n^ . Their magnitude can be estimated by comparing them wltii 
the gravity gradient forces of equation (lJ-28) . Tliese are seen to be propor- 
tional to the square of the term wq . This term, however, is very neai’ly equal 
to the mean angular velocity n. As a consequence, the ineitial forces ai’lsing 
from tlie orbital angular acceleration are by a factor c smaller than tlie gravity 
gradient forces. Because of the extreme smallness of the eccentricity, they 
can be rightfully neglected and were not included in the equations of motion. 
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5. OBLATENESS OF EARTH 


The effects of the oblateness of the Earth on the dynamic response of a 
tethered subsatellite are threefold; two aspects are geometrical in nature, 
while the third is dynamical. The first geometrical effect, the flattening of the 
Earth at its poles, brings about variations in orbital altitude of approximately 
20 km. It has its greatest effect on satellites in polar orbits. The second 
geometric effect results because the orbits of the main satellite are no longer 
circular. In fact, they are not even closed but process in a westerly direction. 
The rate of precession v/ill depend upon the inclination of the orbit and, to a 
somewhat smaller extent, on its altitude. HoAvever, the orbit regression amounts 
to only a few degrees per day and will not induce any discernible dynamic effects 
on the tethered subsatellite. The altitude variations resulting from these non- 
circular orbits range from 5 to 10 Itm, depending on the inclination of the orbit. 

Like the geometric effect of the orbital eccentricity, the geometric 
effects of the Earth’s oblateness wiU result in a constantly changing aero- 
dynamic disturbance force. A detailed simulation of this condition was not 
attempted. It is expected tliat the influence of these effects on the dynamic 
behavior of the tethered subsatellite is very similar to that observed for the 
ox’bital eccentricity. 

The dynamical effect of the oblateness of the Earth is due to the deviation 
of the gravitational field from its spherical symmetry. To determine this effect, 
we begin with the oblateness force exerted on a particle of mass m . Resolved 
in components relative to the orbital reference frame, they are 


Rjj 

2m J g\ ) sin^i sin ji cos p 


= 2m J sin i cos i sin p 


/n \ 4 

/ \ / 1 \ 
3m J g'l ~ ) (3 - sin^i sin p\ 


(11-44) 
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where the oblateness term J = 1,637 x 10"^. The angle /? defines the position 
of the main satellite relative to the intersection of the orbital and equatorial 
planes. The position of any mass element dm of the tethered subsatellite in 
terms of (Ro, can be approximately expressed as 


A) 


and 


/3 + 


R 


1 ^ J_ 

rJ R^ 



(11-45) 


(II-46) 


Analogous to the gravity forces, only the gradient effects of the oblateness 
forces are of interest. These are obtained by inserting equation (11-45) and 
equation (U-46) into equation (11-44) and subtracting the xeroth oblatencss 
contributions. Retaining only linear terms then yields the first-order oblateness 
gradient forces as 


dF 

^0 


R, 


I 


2 2 /3 + 2z sin^ /?) dm 


dF 

^0 


dF = 3 J 


2 wj jI sin i cos i (x cos />’ + 4z sin fS) dm 


"0 


r A 9 9 \ 

\~R~ / r V3 " ^ / 


0 \ R 


siir i sin'^ /i) + x sin^ i sin'^ p 


(11-47) 


dm 


A comparison of these terms with the corresponding gravity gradient terms of 
equation ()I-3) will readily reveal that they are sevoi’al orders of magnitude 
smaller, Therefore, the dynamical effect of the oblatenoss of the Fartii can bg 
s afely di s regarded. 
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6. VISCOELASTIC TETHER FORCE 


The viscoelastic properties of the tether were modeled by very idealized 
elements. Accordingly, the elasticity was represented by a linear spring whose 
spring constant k is 



where d is the;diameter of the tether and E is Young' s modulus. Likewise, 
the damping, whi'eji is a I'ather complex phenomenon, was represented by an 
equivalent viscous damper with a damping foi’ce proportional to the velocity of 
the generalized modal coordinate (modal damp’ug). Consequently, the general- 
ized force due to the elasticity and damping of the tether can be directly obtained 
as 


= - k ^ - Cj C (11-49) 

where C 2 is the coefficient of viscous damping. It is important te h&te that the 
energy dissipation due to frictional losses in the tether material is, in general, 
too small to aid in damping out transient responses of the tethered subsatellite. 
As la pointed out in other sections of this repoi’t the removal of these transients 
requires the employment of an active teth.. — tiol law. 


H. Tether Control Laws 

The dfeployment, retrieval, and quite probably the stationlieeping phase 
of a tethered subsatellite require an active tetlier control system to guarantee 
adequate dynamic \systems performance. In accordance with the principles of 
control system desig'n,' the tether control law will Incorporate information about 
the dynamic state of the system in a closed-loop mode. This can be implemented 
in various ways depending on whether the design is based upon classical or 
optimal control theory. In subsequent chapters, a tether contix)! law is analyzed 
wliich generates a tether tension, T, through a linear feedback law of the form; 


T = kj f + Cj i r kj f ^ 


(11-50) 
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in which the gain values kp Cj, and kg are propei'ly chosen. Implementation of 
this contral law requires the measurement of the tether tension, the rate of 
tether deployment, and the length of the tether itself. The present discussion 
examines the possibility of employing a control law that does not require the 
measurement of the tether tension. Besides being simple, this control law is 
useful for exposing some distinctive features of the dynamic behavior of the 
tethered subsatellite during' its deployment. In this discussion, some facts will 
also be presented on the problems associated with the retrieval of a tetlieTed 
subsatellite. The tether control law to be analyzed consists of the following 
three phasjggt^,^ 

1. ^’Exponential Increase of Deployment Rate 


f = » H for fo < f fi (IT-51 a) 

2. Constant Deployment Rate 

i = c for fi < f fg (II-51b) 


3. Ex|X)nontial Decrease of Deployment Rate 


i = 0'(f^ I- fo - n for fg < ^ 

Numerical values used for the dynamic analysis are; 


to " 10 m 

tg = 9 X K)-' m 

c = 5 m/s 


fi ^ 10'* m 

f = 9. 999 X lO"* m 
c 

O' - 5x s"' 


(ll-51c) 


Time histories of the tether length and its deployment rates based uj^on this 
so-called exponential central law are given in Figures II-3 and 11—1. To gain 
physical insight into the dynamic behavior of the system, the attitude equations 
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were examined for both the in-plane and out-of-plane motion in the 
absence of aerodynamic disturbances. Furthermore, the tether mass was 
assumed to be negligible and the orbit of the main satellite to be circular. For 
this special' case, the in-plane motion is decoupled from tlie out-of-plane motion 
and is described by 


0 + 




3 ? . „ 

+ — Wq sm 2 0 



(11-52) 


The out-of-plane motion only decouples from the in-plane motion for small pitch 
angles 0 . With this assumption, the out-of-plane motion is described by 


(p + 2' 


(f)^ 


4 Wj 0 =0 


(11-53) 


At first tlie dynamic behartor of the system is examitied for constant tether 
length (I = 0) . For this case, both the in-plane and out-of-plane motion are 
analogous to a simple pendulum motion with stable equilibria at 0 = O’ and 
0 = 180° . However, the in-plane oscillation has a different period than the 
out-of-plane oscillation, even for small attitude angles. They are given as 

Out-of-Plane Period 


To = 


IT 

n'o 


( 11 - 54 ) 


In-Plane Period 


2 IT 
Wo 


1.15 To 


The in-plane period is seen to be somewliat larger than the out-of-plane period. 
It is important to nocice that the periods of both oscillations are independent of 
the tether length. 
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For large attitude angles, the in-plane and out-of-plane motions are no 
longer decoupled. As a consequence, the tethei’ed subsatellite will, in general, 
exhibit a rather complex motion. Because the frequencies of tlie two oscilla- 
tions are incommensurable, the resulting pattern of the motion will somewhat 
resemble the so-called Lissajous figui’es encountered in celestial mechanics. 

During the deployment and retrieval maneuvers, the tether length is no 
longer constant. Examination of the preceding equations of motion reveals that 
13ie dynamic behavior of the tethered subsatellite becomes castinctively different 
for these two maneuvers. For deployment, the tether lengh I is increasing 
(f > 0) and the second tei-ms of equations (11-52) and (11-53) containing the 
derivative of the tether length tal^e the place of an equivalent damping. In such 
a case, existing oscillations will have a tendency to die. I'or retrieval, the 
tether length is decreasing (t < 0) and the second term introduces negative 
damping. In this case, existing oscillations will tend to build up. It has been 
fomid b}' simulation that this pattern is very pronounced in the out-of-plane 
motion. Because of the presence of the term of the right side of equation (II-52) , 
however, the in-plane motion does not build up rapidly. It is expected that this 
accumulation of kinetic energy during tlie retrieval phase can be held within 
acceptable limits by a proper design of the tether control lawn Such a control 
law will, by necessity, lead to a "yo-yo” type of retrieval motion. This motion 
wiU be characterized by fast retrieval rates in the neighborhood of tlie maximum 
attitude excursions of the subsatellite and slow or even nt gative retrieval rates 
(i. e. , i >0) in the vicinity of the equilibrium points. Bi. causc of this yo-yo 
effect, the retrieval time of a tetliered subsatellite will, oeteiis paribus , be 
longer than tlie time for its deployment. Some preliminary studies concentrating 
on the problems associated with the retrieval of a tetherr d subsatellite are pre- 
sented in a subsequent section. One can conceivably als': use this yo-yo effect 
for the removal of residual kinetic energy finm a tetherc d subsatellite during 
its stationkeeping' mode. An optimal bang-b;mg type, yo -yo control law for this 
purpose is illustrated in Figure n-5. Each time the sub satellite passes through 
its equilibrium position the tether length is suddenly decreased by a small 
amount 6 and then increased by the same amount when ;he subsatellite reaches 
its extreme position. The dissipation of Idnetic energy ))er cycle resulting from 
this maneuver can be related to an equivalent viscous damping effect. For small 
oscillations of the in-plane motion, the corresponding c itical damping ratio 
becomes 




2 / 6 \ 

7T f ) 


(11-55) 
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where Oo is the amplitude of the in -plane oscillation. It is seen that the damp- 
ing effect increases with decreasing amplitude, indicating that yo-yo control of 
the In-plane motion could be very effective. The corresponding damping effect 
for the out-of-plane motion can be calculated as 





(11-56) 


This damping effect is independent of the amplitude, indicating that a yo-yo 
control of the out-of-plane motion could be very ineffective. In fact, a yo-yo 
motion amounting to a 3 percent change in tether length induces only a 1 percent 
damping ratio. 

A schematic presentation of this yo-yo control law is shown in Figure 
II-5. It should be realized that the bang-bang type operation of this optimal 
yo-yo control law prevents its direct hardware implementation. It was only 
discussed for demonstrating the dynamic principle involved in removing kinetic 
energy from a system without the expulsion of mass by reaction jets. 

Inspection of the in-plane motion, ecmation (n-52) , also uncovered an 
interesting and important dynamic peculiarity associated with the deployment 
and retrieval phase of the tethered subsatellite. It is seen that the subsatellite 
can assume a steady-state attitude angle 0 if tlie deployment velocity increases 
exponentially as in equation (ll-51a). However, stability of the motion can only 
be maintained as long as the attitude angle stays below 45°. Obviously, this 
limits the deployment rate a , which, in turn, establishes a lower bound for 
the deployment time of the subsatellite. For negligible tether mass, the upper 
limit of the deployment rate, , is constant with time. If the tether mass is 

taken into account, the system reaches only a quasi-steady state and the upper 
limit of the deployment rate changes witli time according to the relationship 
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where m=m(t) is the instantaneous tether mass. To verify this result, com- 
puter simulations were made with the deployment rate a as a parameter. The 
dynamic behavior of the subsatellite during the exponential buildup phase of the 
deployment velocity is sho\vn in Figure II-6. The critical deployment rate was 
determined by iteration to be approximately « 9. 5 x 10"^ s“*. This value 

agrees rather well with the quasi-steady value of equation (II-56). Physically, 
the onset of the instability is caused by the loss of tension in the tether. From 
this time point on, the subsatellite becomes a free-flying object imtU the ■*:ension 
is regained. For all practical purposes, however, loss of tension in the tether 
also signals loss of control over the subsatellite. The same dynamic peculiarity 
also exists, of course, for the retrieval phase with the proper sign I’eversal of 
tlie attitude angle 0 . Instead of trying to stay "in front" of the main satellite as 
during its deployment, the tetliered subsatellite will now stay "behind" the main 
satellite (o < 0) . As was discussed preriously, the retrieval phase is, however, 
dominated by the second term of equations (11-52) and (11-53) , especially for 
the out-of-plane motion. The negative damping effect caused by tliis term 
becomes very noticeable in the vicinity of tlie main satellite where the tether 
length is small. There the tethered subsateUite acquires a spectacular tendency 
to wrap around the main satellite, imless the retrieval rate is reduced to an 
extremely small level. It should be pointed out that these dynamic characteris- 
tics of a tethered subsatellite during its deployment and retrieval are intrinsic 
in nature and therefore independent of the particular tether control law chosen 
for these maneuvers. 


1. Results of Computer Simulations 

To gain a conclusive engineering imdei standing of the dynamic behavior 
of the tethered subsatellite system, the analytical investigations were supple- 
mented by detailed computer simulations. These were based on design specifi- 
cations proposed by die Smithsonian Astropliysical Observatoiy. Accordingly, 
the following system parameters were used: 

Main Satellite 


Orbital Altitude 
SubsateUite 

Spherical lUass 
Cross Section 
Drag Coefficient 


n 200 ion 

M - 170 kg 
A - 1 m^ 
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Tether 


Deployed LengUi 

t = 100 Ian 

Radius 

r = 0.183X 10“5 m 

Density 

p = 7. 93 X 10^ leg m"^ 

Total Mass 

m = 83. 4 kg 

Normal Force 
Coefficient 


Axial Force Coefficient 

C. = 0.2 
A 


The sensitivity of the dynamic response characteristics was investigated with 
respect to variations of both system parameters and initial conditions. This 
revealed a number of interesting and important aspects of the system. First, 
it was noticed that the system is relatively insensitive with respect to variations 
of the initial deployment conditions. Dynamic transients ai’ising from these do 
not build up but stay within acceptable limits. The initial alignment of the sub- 
satellite and its actual release mechanism should, therefore, not be a critical 
item. Likewise, it was foimd that the dynamic response characteristics were 
not very sensitive to changes in payload and tether mass. The design of the 
tethered subsatellite system can, therefore, accommodate a rather wide range 
of these parameters. On the other hand, the effect of orbital eccentricity was 
very pronomiced because of the drastic atmospheric density variations induced 
by the altitude changes of the subsatellite. Care has to be taken, therefore, to 
maintain a nearly circular orbit of the main satellite during the operation of the 
tethered subsatellite. For a subsatellite altitude of 100 km, it seems advisable 
to limit tlie altitude variations due to orbital eccentricity to 10 km or less. 

This requirement becomes, of course, less stringent for higher subsatellite 
altitudes because of the decreasing effect of the atmospheric disturbances. 

To illustrate the d 3 mamic behavior of the tetliered subsatellite, only a 
few typical cases wei’e selected, and these are shown in Figures n-7 through 
n-20. The plots are arranged in the following order; 

1. Pitch Angle 0 

2. Roll Angle 0 
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3. In-Plane Motion (x - z Plot) 

4. Out-of-Plane Motion (y-z Plot) 

5. Equai rial Motion (x-y Plot) 

6. Tether Elongation ^ 

7. Tether Tension T 

Initial conditions for all cases are 0 q = %= 3“ and = 0 , = wq . They 

reflect initial misalignments and a small out-of-plane release impulse of the 
subsatellite. Each figure contains two plots: one for a circular orbit (Case A) 
and one for an eccentric orbit with AH = 10 km (Case B) . Figures 11-7 through 
11-13 present equatorial orbits and Figures 11-14 through 11-20 present polar 
orbits. Cases of intermediate orbital inclinations were not included because the 
general nature of their dynamic response characteristics can be easily visualized 
by proper inteipolation between the cases shown. The time histories of the state 
variables presented are essentially self-explanatory. However, a few remarks 
seem to be in order. For equatorial orbits, the out-of-plane motion (y-z plot) 
is essentially dictated by tim initial conditions. Since the exponential deployment 
does not provide damping, the responses are very nearly constant. The equa- 
torial motion (x-y plot) shows the typical Lissajou characteristics, especially 
for the fully deployed condition. Eccentricity of the orbit introduces substantial 
excursions of the in-plane motion (x-z plot) with dynamic transients of several 
thousand meters over and above the steady-state response. For polar orbits 
the effect of the rotating atmosphere is clearly visible. It leads to out-of-plane 
transients of 2000 to 5000 m. However, there are no indications of resonance. 
The tensile forces in the tether never exceed 100 N. They should not present 
problems relative to tlie material properties of the tether. The same holds true 
for the elastic expansion of the tether w'hich stays below 400 m or 0, 4 percent 
of the total tether length. 


J. Steady-State Lateral Tether Deflection 

Because of limitations in manpower and resources, only a cursory 
examination could be made of the effect of lateral tetlier dynamics. The follow- 
ing analysis makes use of the Rayleigh-Ritz method, which is based on the 
premise that the exact deflection of an elastic system can be approximated by 
the superposition of suitably chosen mode functions. The success of the 
Rayleigh-Ritz 2 netliod depends very much on the choice of these mode functions. 
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They are usually selected from tlie natural mode shapes (eigenfunctions) of the 
dynamical system under investigation or of one which is very similar to it. For 
the tetliered subsatellite system, it wo\dd be natural to use the eigenfunction of 
a string vibrating under nommiform tension. This would be determined by the 
gravity gradient force field. The string would be suspended at one end and have 
a mass attached to the other end. Since a computer program was not readily 
available to compute these eigenfunctions, the analysis was based on eigen- 
functions of a string under uniform tension. The string was assumed to have 
an infinite end mass. To provide for maximum tetlier flexibility, the tension 
was set equal to the smallest value of the tension in the tether. This occurs at 
the attach point of the subsatellite. For this condition, the eigenfunctions, tp , 
are simply given by 


t! (z) = sin 
n 


n 7T z 


n = 1, 2, . . . 


( 11 - 58 ) 


The associated eigenvalues (natural frequencies) of tlK system are 


CO = 


n 



( 11 - 59 ) 


In accordance with the Rayleigh-Ritz method, the modal response of the system 
is obtained from the equation 


n + 
'n 


2 

W ^ 71 

n 'r 




M 

n 


( 11 - 60 ) 


where f(4>t) represents the aerodynamic force per unit length of the tether 
7 )^ the generalized coordinate, and the generalized mass of the nth mode 

shape. The total response of the tether over and above the rigid tether I'esponse 
is then obtained by supei’position of the individual modal responses as 


u(z,t) =-• T i/'j(z) J)j(t) i = 1, 2, ... N 


( 11 - 61 ) 
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where N is the total number of mode shapes used in the Rayleigh-Ritz approxi- 
mation. For the present feasibility study, the response of the fully deployed 
tether was only detennined for a steady-state aerodynamic excitation force. 

This very closely corresponds to a tethered subsatellite moving in a circular 
equatorial orbit. For eccentric or inclined orbits, the aerodynamic forcing 
function contains a component that has an approximate sinusoidal variation with 
time and a period close to the orbital period of the orbiter. However, it should 
be remarked that the period of the first few tether mode shapes is vei’y much 
smaller than the orbital period. Therefore, no serious resonance phenomena 
are expected to occur for these orbits, at least not for the fully deployed con- 
dition. 


The total number N that has to be used to obtain a satisfactory approxi- 
mation to the actual deflection of a dj namic system depends strongly on the 
forcing function itself. Because of the drastic exponential Increase of the 
atmospheric density with decreasing altitude, the aerodynamic excitation is 
concentrated at the lower end of the tettier. As a consequence of this strong 
nonuntfomi force distribution, it takes a relatively large number of mode shapes 
for convergence. This fact is illustrated in Figures 11-21 through II-25. It is 
seen that satisfactory convergence requires the inclusion of 8 to 10 mode shapes. 
The maximum deflection of the tether is about G km, i. e. , about 6 percent of 
its total length. 

A more accurate analysis of the lateral tether flexibility would have to 
be based on Lagrange's equation (11-9) using the previously mentioned refined 
mode shapes. 


K. Conclusions and Recommendations 

The feasibility study conducted thus far did not expose any hazardous or 
adverse dynamic response cliaracteristics of the tethered subsatellitc during its 
deployment and statlonlceeping j^hase. In fact, its dynamic behavior was imex- 
pectedly favorable. Part of this can, of course, be attributed to the fact that a 
rather long duration was allowed for the deployment maneuver whose initial 
velocities wore kept at a few millimeters per second. In practice, tliis would 
lead to deployment times of 8 to 12 h. However, there seemed to be no objection 
to allotting this amount of time to the total deployment pliase. The dynmnic 
transients and steady-state excursions during the stationlvccping phase appeared 
to be acceptable for the proposed satcllite-bome experiments. Besides, there 
seems to be tlic possil^ility of reducing these dynamic responses by a properly 
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chosen yo-yo tether control law. This could ob\iate the potential necessity of 
having an active reaction-jet control system on board the subsatellite. Tensile 
forces in the tether were well within the strengih limits of the proposed tether 
materials and gauges. No dynamic resonance phenomena were observed as a 
result of aerodynamic disturbances for inclined or eccentric orbits. Although 
diurnal fluctuations of the atmospheric density were not simulated, they will 
probably not introduce serious dynamic response problems. They will, liowever, 
slightly increase the steady-state excursions of the subsatellite. Whei’eas the 
dynamical effects of the Earth's oblateness were found to be negligible, its geo- 
metric effects are significant. A satisfactory assessment of these would have 
to be based on a more accurate definition of the atmospheric density variations 
above the surface of the Eai-th including diurnal effects. It did not come as a 
surprise that a serious dynamic problem arose during the attempt to retrieve 
the subsatellite. This problem is characterized by a spectacular wrap-around 
phenomenon during the last phase of the retrieval maneuver, when the sub- 
satellites come into tire close vicinity of tire main satellite. A similar behavior 
had already been observed in earlier dynamic studies that were performed to 
prove the feasibility of rescuing a disabled astronaut by tether retrieval. 

Since this problem and its potential solution will be discussed in a later section 
of tills report, it will not be pursued furiher here. It seems appropriate to 
mention, however, that an earlier planar analysis of the retrieval phase did not 
reveal this phenomenon because the tether control law could provide sufficient 
damping for the in-plane motion. Therefore, it is recommended that future 
analyses be performed using a three-dimensional dynamical model of the 
tethered subsatellite. 

Based on the physical insight gained in the present dynamic analysis, it 
is furHier recommended that the dynamic analysis be expanded by including 
lateral tether mode shapes in both pitch plane and yaw plane. In addition, tiie 
effect of the tethered subsatellite on the main satellite should be assessed, in 
order to estimate the propellant consumption for maintaining proper attitude 
and orbital position of the main satellite. This can be done in a separate analysis. 
Some effort should be devoted to analyze the effect of the ionospheric plasma and 
the Earth's magnetic field on the tether. Both of these effects are expected to 
be small. 
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Figure II-3. Exponential deployment law — rate of tether deployment versus time. 
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Figure II-ll. Equatorial motion (x-y plot) — equatorial orbit. 
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Figure 11-13. Tether tension versus time — equatorial orbit. 
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Figure 11-14. Pitch angle vex’sus time — polar orbit, 
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Figure 11-19. Tether elongation versus time — polar orbit. 
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Figure n-21. Steady state modal deflections — 1 and 2 modes. 






Figiare 11-22. Steady state modal deflections — 3 and 4 modes. 
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Fig-ure H-23. Steady state modal deflections — 5 and 6 modes. 
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Figure 11-24. Steady state modal deflections — 7 and 8 modes. 


I 

cn 




DEFLECTION (km) 



Figure 11-25. Steady state modai deflections — 9 and 10 modes. 




III. INVESTIGATION OF A TETHERED 
SUB SATELLITE CONTROL LAW 


Ralph R. Kissel 

A. Summary 

A control law is presented which uses tension in the line, line velocity, 
and length to deploy and retrieve a subsatellite up to 100 Ion above or below 
the Shuttle or other space vehicle. The equations of motion used here are 
three-dimensional but include only the mass of the line and no other factors, 
i. e. , tetlier dynamics, atmospheric forces, elliptical orbits, etc. These addi- 
tional factors are covered in detail by M. Rheinfurth in Section II. 

The results show that tlie subsatellite can be deployed in a stable manner 
in a few hours and can be reti’ieved in a stable maimer if enough time is avail- 
able. Approximate^ 19 h is the minimum time required for retrieval to 5 m 
from 100 1cm. 

The subsatellite can be positioned to any distance from any other dis- 
tance and stabilized at the new position in a few hours (107 up to 120 ±0.5 Icm 
altitude required 1 . 5 h) . 

Feasibility of tliis technique will be governed by mechanical design 
limitations and time allowed to stabilize at a new position. 


R. Introduction 

Until the work presented in this report was done, no continuous control 
law was available to deploy, retrieve, or position a subsatcllite suspended by a 
tether up to 100 tan long. This law is an extension of one developed by C. Rupp 
and used by AMA in their work. The extension involves maldng the commanded 
length in that original control law a function of the actual length rather than an 
arbitrary set of commands. The particular function chosen is determined by the 
desired line velocity, position, damping, and retrieval time. 


m-1 



C. Method 


The dynamic equations used here were derived by M. Rheinfurth and 
are presented elsewhere in this publication. The control law developed by 
Rupp is 


T = 7a,'o^i + 4o)oi - 


(m-1) 


where 

0)0 = orbital rate 

t = tether length 

S = is commanded length 
c 

T - commanded tether accelei’ation. 


The tension in the tether ( at Shuttle) is given by 


T = 



M 

wire 


(m- 2 ) 


where M is the subsatellito mass and M . is tlie mass of wire of Icngtli f . 
sat wire 

For purposes of optimizing the tension, the control law was expanded to 

T = (R^ + 3) + Zijiwof - (ni-3) 

where R is the oj /oin. w is the control law "stretch" frequency, and f is 
c c c 

the control law damping. 


m -2 



The control law extension developed here sets 


= Kii + K 3 (m-4) 

where Kj is chosen for performance and K 3 is chosen for the desired deployment 
velocity. 

The idea for this extension was based on the old "carrot-and-stick" idea, 
i. e. , command a length always ahead of the present length by an amount which 
will give the desired trajectory. 

The values used for Kj and K 3 are different for retrieval and deployment. 
During deployment, three sets of values are used depending on the relationship 
of the current length to either the phase 1-2 mode change length or the final 
length. During the initial phase, Kj is given a value greater than one while K 3 
is set to zero. The value of Kj for this phase is chosen to minimize the time 
to the desired velocity. Wlien the desired velocity is reached, the second phase 
begins and continues until the terminal phase begins. During die second phase, 
is set to one and K 3 is chosen to maintain the desired velocity} A smooth 
transition from the initial phase to the second phase requires that 


K 3 (phase 2) = L[ Kj( phase 1) - 1] 


(m-5) 


where L is the length for changing from phase one to phase two. The terminal 
phase has K( = 0 and K 3 equal to the desired final length. The terminal phase 
begins when the actual length is within K 3 (phase 2 ) units of the final length. 

Retrieval works best when K 3 = 0 and Kj < 1 for the entire retrieval time. 
If Ki is too low, the subsatellite will go liy the Sliuttle; if K] ajiproachcs one then 
too much time is required to retrieve. 

Other parameters which can be optimized together with K] and K 3 are 

(u and t . The w and f ai’c chosen primarily to minimize out-of-plane 
c c c c 

motion during retrieval. /Vn cu of four coo is theoretically optimum for best 
damping of out-of-plane motion. 
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station keeping is dependent on w and ^ , but not on Kj. The values 

c c 

for 03^ and ^ ^ will need to be chosen based on desired response and may be 

different from values used during retrieval. The value for Kj is similar to 
that used during the third phase of deplo 5 onent, i. e. , it is the desired length. 

Implementing this control law may best be accomplished by a combina- 
tion of digital and analog teclmiques. Equation (ni-l) can probably best be 
implemented using analog teclmiques while the choosing of Kj and K 3 appears 
simplest by digital metliods. 

Mechanizing the entire system has several difficulties, most of which 
must still be fully resolved. These will be briefly disrusscd in the conclusions. 

All the simulations for this report were done on a desk-top HP-9830 
calculator. Ail programs were either written especially for this problem or 
were adapted from an earlier project involving similar requirements. 


D. Results 

Figures III-l through ni-6 are illustrations of what could be expected. 
Obviously, a large number of simulations could be done depending on the control 
desired and the mission to be accomplished. Once those are specified a partic- 
ular optimum solution could be worked out. 

Table ni-1 shows the parameter values used for the simulations in this 
report. Also used wci’o 


M , = 170 leg 
sat 

M . = 83 kg (100 1cm) 

wire 

= 2 Wo 




c 


1 


with coordinates as defined in the list of symbols. 
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TABLE m-1. INITIAL CONDITIONS AND CONTROL VALUES 


Units 

m 

m/ s 

deg 

deg/ s 

deg 

deg-/ s 

■None 

m 

Name 


^ 0 

So 

O'o 

i/'O 


Ki 

Ks 

Deploy 

10 

0. 5 

0 

0 

0 

0 

2 

7 000 

Retrieve 

100 000 

0 

20 

0 

10 

0 

0.93 

0 

Station Keeping 

100 000 


0 

20 

0 

10 



0 

0 

10C5 000 




























Figure III-l shows a side view of the deployment and retrieval trajec- 
tories. Deployment is smooth and relatively fast, with average tether velocity 
being purposely limited to approximately 8 m/ s. Tliis limit was deemed con- 
sistent with reel capability in steady state. The resultant time was 6. 7 h to 
stabilize at f = 100 Ian (atmosphere not considered) . Setting the subsatellite 
out to 20 m with no initial velocity I’equires only 15 min more to go to 100 Ian. 
There was a 3 Ian overshoot but a change of parameters during the terminal 
phase could reduce this to almost any desired value. 

Retrr’val began with an extreme worst case, i. e. , with the subsatellite 
in front of the Shuttle path and slightly off to the side. Atmosphere would tend 
to pull the subsatellite behind the path which other simulations show does not 
cause such a low altitude to be reached as shown here. Different parameters 
can minimize this effect to almost any point desired although somewhat more 
time would be required and pai-ameters would need to be changed during retrieval. 
This particular retrieval required 26 h to reach 5 m below the Shuttle. The 
effect of out-of-planc motion is seen towards the end of tlie retrieval by the 
slight in-plane component. 

Figure III-2 shows the same retrieval from the front. The out~pf-planc 
motion is clearly seen. It can also be seen how this motion is dami^ing put as 
the subsatollite approaches the Shuttle. 

Figure III-3 shows the tether tension and velocity throughout deployment; 
both are relatively smooth. Maximum tension at the Shuttle is approximately 
93 N and maximum velocity is approximately 9 m/ s. 

Figure III-4 shows tether tension and velocity during retrieval. Maxi- 
mum tension at the Shuttle for this extreme case is about 123 N and maximum 
velocity is approximately 23 m/ s. Velocities vary (plus and minus) and tension 
and velocity oscillate as a result of the control law tri'ing to dampen the motion. 

Table III-2 shows the best retrieval times to a given distance below the 
Shuttle that can be oi 'alned with this control law. 


TABI.E m-2. BEST RETRIEVAL 'ITMES TO f 


i (m) 

5 

10 

50 

100 

t(h) 

19,1 

18.8 

17.2 

10. 1 
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These times can be obtained with any I and o) but their choice is 

c c 

limited by the oyejj^hoot and maximum line velocity allowed. Times shorter 
than those in Tjtbl© 111-2 cause the subsatellite to miss the Sliuttle. 

Thb JaSt part of retrieval requires precise iaiowledge of the control law 
parameters £, £ , and t. Table IIfr-3 show^- 1 ^, resolutions required for auto- 
mated retrieval to 50 ni. 


TABliE in-3. RESOLUTIONS FOR AUTOMATED 
RETRIEVAL TO 50 m 


Parameter 

Range 

Resolution 

Length 

50 to 102000 m 

±5 in 

Velocity 

-16 to 22 m/ s 

viO. 005 m/ s 

Tension 

0.001 to 130 N 

i0.0025 N 


Station keeping is illustrated by two very different cases in Figures III-5 
and in-6. Figure III-5 uses the parameters given in Table rU-l and on page 
in-4, and Figure UI-6 shows results by changing;:j*f^v >^ 4 wq and = 

10. The standard parameters allow the tether length to a)5|)roach 114 km and 
after 2 h is within 0. 5 km of the desired 100 Ion length. This case clearly 
could not be tolerated. Figure III-G shows the opposite situation in that the 
line length never e: - eded 101 Ian, and since this length did not occur directly 
beneath the Shuttle, the subsatellitc never even exceeded 100 Ian below the 
Shuttle. The actual parameter’s for station keeping would necessarily be 
between these two extremes. 

Best station keeping values (for best daminng) arc w = 4 wq and ^ = 

0. 5. Approximatclj' 8 h arc needed to cut the out-of-plane amplitude in half 
after an initial disturbance. Velocity and tension are similar to that cxpei'i- 
enced during the initia,! part of retrieval . 
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FigT^re III-7 shows a possible analog diagram of the reel control system . 
The control law produces^ttg desired torque and this is compared with the torque 
actually existing, measui’|^n this case at the motor, to produce the error 
signal. 


Figure in-8 shows the digital logic needed to implemei^^thc selection of 
Kj and K 3 for each phase of the trajectox'y. The desired lengtiws compared 
with the last commanded length and the current length to make the pi'oper 
selection. Also needed is the predetermined value for changing from piiase 
one to phase tv'o in deployment. 

Figure III -9 shows more detail on an actual mechani’S^^B'on. The torque 
motor uses tachometer feedback for stabilization. The length and velocity of 
the tether are measured by a pulley kept in frictional contact vith the tether. 
The tension is measured by a spring damper arrangement on this same pulley. 
These signals are fed to the processor to produce a torque command. 

Other ideas have been put forward as possible ways to deploy, retrieve, 
and control the subsatellite. Initial deployment methods include a spring 
release system, a gas jet driven system, and a rail or tube guided system. 
Retrieval ideas include using a boom and hook, a cushioned cone or catcher, 
or moving the Shuttle domi to the subsatelHtc. 

Length measurement could bo done using magnetic jDulscs on the wire, 
paint, and photo-optics readout, counting turns of the wire reel or counting 
revolutions of a pulley hold against the wire. An eccentric motor gear would 
give continuous measurement. 

Direct position mcasm-ements may be possible by radar, lasci , strobe 
light, received power, or measuring tether angle and length relative to tne 
Shuttle. Some of these arc discussed in other parts of this publication. 

Wire velocity moasurctnont could use the length measurement ideas and 
get a differential measurement per unit time. 

Tcnsi'caj/ieould be measured by the spring system already mentioned or 
by ineasuring motor current. 
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The required motor power has been estimated at approximately 900 W 
peak to develop about 9 N m at up to 600 rpm. Such a motor weighs approxi- 
mately 50 N (10 lb) . For a full 0.6 m diameter reel weigliing about 5800 N 
(1300 lb) this motor could produce about 0. 12 m/ s^ linear acceleration. This 
is rdequate for all control requirements. It would also be adequate for initial 
deployment if 4 or 5 s is allowed for reel spin up or if a teclmique is used wMch 
does not require Mgh reel acceleration. 


E. Conclusions 

Deployment and i-etrieval are theoretically feasible based on these con- 
trol law studies. There are several potential problem areas, however. 

First, the long x’etricval time ma^' limit the usefulne.ss of the concept. 
Some time can be saved by retrieving' the subsaf^ellite at some distance from 
the Shuttle if a m-'uiod is worked out to do tills, such as a long' boom or moving 
the Shuttle. 

Second, retrieval measurements just before retrieval occurs require 
high resolution to guarantee success. A single device to measure tension, 
length, and velocitj' to the required I’esolutior. over their resjiective ranges 
could be difficult to build. However, a scheme to diidde the range into two or 
more parts could reduce this difficulty. 

Th.- mechanics of actually deplojdng and retrieving the subsatellite 
could be quite involved. The subsatellite must be either gi'/on an initial velocity 
or displaced from the Shutllc for gravitj' gradient to become effective for further 
deplojonent and control; several ideas have been advanced to do this. Retrieval 
meehanisms must allow for a certain error in both velocity and position at 
retrieval time. They must also retrieve without damage and with guaranteed 
success; ideas have also been presented to do tliis. 

Methods of mea.suiing length, velocitj’, and tension have been advanced, 
but the merits of each must still be evaluated. 

High velocities w'ere usually seen at some time during retrieval, espec- 
ially if minimum retrieval time was being sought. This places added require- 
ments on the reel co’'.r.'>t mechanism ant cable construcnoi . 
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The trade-off for the best control law for a given mission must be done 
whenever the mission is specified. This may require, in the most complicated 
case, that all control parameters be functions of time. Station keeping, 
retrieval, and deployment would each have their own parameter set. 

Operation above the Shuttle is the same as operation below'. Missions 
being mentioned thus far require a different length and size cable for operation 
above than below'. Tliis could mean another control parameter set, at worst, 
but at least this w'ould have to be investigated further. 

Maximum tether acceleration has been less than 0. 1 m/s^ w'liich is not 
a problem. How'ever, initial deployment, if an initial velocity is required, 
could produce a far liigher value. It is this high value that sizes the reel con- 
trol motor if the entire reel must be accelerated; tliis must be furtlier worked 
out. 


The control law' presented here allows freedom to stop and start from 
any tether length. This simply says that the control law’ is continuous rather 
than stepwise or arbitrary. Tension measurement is a requirement for suc- 
cessful retrieval. It is not a necessity for deployment. 

Further work should be done to test this contro, law' in the detailed 
simulation explained elsewhere in tliis publication. Also a liardware simulation 
could be of considerable value once a specific mission is defined. 
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Figure IH-2, Retrieval front view. 
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Figure III-5. Station keeping front view. 
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Figure m-8. Flow chart to find Kj, Kj for ^ Ki* i + K3. 
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Figure III-9. Conceptual mechanization. 


iV, AERODYfJAMlCS 

Kenneth D. Johnston 


The aerodynamic forces and heating' on the satellite and tether wire are 
estimated and one possible design for aerodynamic attitude stabilization is 
examined in this section. The interference of the satellite -tthth atmospheric 
probe measurements is also briefly discussed. 

A. Tether Wire 

Tlie tether wire is so small (D ^ 0.4 mm) that the wire is in the free 
molecular flow regime over the entire altitude range (altitude 5^ 100 km) . The 
Knudsen number (molecular mean free path/ wire diameter) of the wire at 
100 km altitude is approximately 400, and it increases at higher altitudes. A 
Knudsen number of 10 or greater is sufficient to guarantee free molecular ilow. 

The aerodynamic force coefficients were estimated using Reference IV-1. 
Although the free molecular force coefficients vary somewhat with wall tem- 
perature and speed ratio, the accuracy obtained by assuming these quantities 
to be constant over the altitude range is considered adequate for this study. 
Therefore, these quantities were taken to be constant at the values they assume 
at the midaltitude of 150 km. The normal force and axial force coefficients 
for the wire are given in Figure IV-1 as a function of the angic-of-attack of the 
wire. 


Experimental data from Reference lV-2 were used to estimate the aero- 
dynamic heat input to the wire. Figure IV-2 gives the total film coefficient, 
li,^, (defined in the figure) , as a function of altitude for free molecule flow at 

I’ight angles to the wire. The equilibrium temperature of the wire was estimated 
by balancing the acrodjmamic and solar heat input against the emitted radiant 
energy assuming that both the ab.soi’btivity and omissivity are equal to 0. 85. All 
other sources of energy input such as Earth radiation and electric current in the 
wire were igiioi’ed. Figure IV-.l gives the equilibrium wire temperature as a 
function of altitude with and without solar heating. 
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B. Nonaerodynamically Stabilized Satellites 


If there is no requirement for aerodynamic attitifdo control, simple 
external shapes such as spheres or cylinders can be used. Only spheres and 
cylinders are considered here. Four sizes ranging from 0. 5 to 2 m in diameter 
were chosen to bracket tlie probable size of the subsatellite. The speed of the 
subsatellite was obtained, whenever needed, by observing that the subsatellite 
and the orbiter travel at the same angular speed about the center of the Earth. 

Figure IV -4 gives drag coefficients of spheres as a function of altitude. 
These curves were obtained from experimental data in Reference IV-3 for 
hypersonic, rarefied flow. A different curve is obtained for each diameter at 
the lower altitudes whore the spheres are in the transitional flow I'cgimc, i.c. , 
that flow regime between the continuum regime at low altiUiC'.'^ .md the free 
molecular regime at high altitudes. Moderate changes in boc y si e produce 
significant changes in force coefficients in the transitional flow rci imo. 

The aerodynamic heat input to a 1 m diameter sphere w is estimated 
using experimental data from Reference IV-2. The total film coefficient, li^^,, 

is given in Figure IV-5. An cnei’gy balance similar to that for the tether wire 
was performed to find tlie equilibrium temperature- of the sphere. This tem- 
perature is presented as a function of altitude in Eiguro IV-O. 

Drag coefficients for cylinders at zero angle-of-attack and Icngth-to- 
diameter ratio of 2 are given in Figure IV-7. Again, the separate curves for 
the different body sizes indicate the transitional flbw regime. A "bridging- 
formula" given in Reference IV-4 and modified in Reference IV-5 was used to 
cstinivite the coefficients in the transitional regime. 'J’ljo .bridging formula 
bridges the gap between the continuum and free molecular regimes by supplying 
a curve shape based on experimental data for similar bodies. The free molec- 
ular coefficients were obtained fi’om Reference IV -1. 

C. Aeroclynamicaliy Stabilized Satellites 

U' aerodynamic attitude stability is desired, some system of vanes must 
be placed on the leeward end of the satellite to keep it pointing into the relative 
wind. The force exerted by the tether wire on the subsatellite could produce 
moments which might counteract the aerodynamic moments. This interference 
could be prevented if the wire were attached at the center of gravity of the sub- 
satellite. This effect might also be minimized by locating the wire attachment 
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point judiciously. The tether wire,..\vould liavo mucli less effect in the yaw plane 
than in the pitch plane, so the tetlpi^Vwire might be ignored in a first approxi- 
mation for yawing motion. 

For a fin surface to be effective in high altitude, low density flight, it 
must be inclined at a large angle to the flow. Tins bCibecause the lift curve 
slope, 9 / 9a; , of a fin surface is very small nea^tgidro angle of-attack and 

becomes much larger at greater angles-of-attack. No attempt was made to Lnd 
an optimum vane configuration since the present study is only to determine 
feasibility. Therefore, a simple coi®guration was chosen for convenience of 
computation. This configuration consists of a lightweight 45° half-angle cone 
frustum attached to a 1 m diameter spherical s.atellite with the apex of the 
cone coinciding with the center of the sphere. To estimate the moment of 
inertia of tins configuration, it was assumed that the supporting structure and 
skin of the cone-frustum afterbody are equivalent in mass distribution to a 

0. 005 in. thick aluminum conical shell. Cone base diameters of 4, 6, and 8 m 
were chosen for the study. Some physlc-Tl chni’acteristics of these bodies are 
listed in Figure IV-8. 

1. STATIC AERODYNAMICS 

Tlie drag coefficient of this configuration at ?.ero^! 5 ,gIe-of-attack is 
given in Figure F/-9 as a function of altitude. These dawwere comiiutcd using 
Reference IV -4 and the modified bridging formula of Reference IV-5. The drag 
coefficients are convojjtlid to drag force in Figure IV-10 and compared witli the 
sphere alone drag to illustrate the penalty in drag which results from attaching 
the afterbody. The drag force is rather large at the lower altitudes. This 
drag imposes a horizontal component of tension in the wire which might limit 
the size of afterbody whicli can be used. It would be useful to measure the ten- 
sion in the wire at th^ satellite to provide an estimate of atmospheric density 
by assuming that the drag coefficient is Icnown. 

A "design altitude" of 107 km was selected for tliis study. The axial 
force tocffidicnt is given as a function of angle-of-attack for this altitude in 
Figure IV-11. The normal force coefficient is given in Figure IV-12. The 
aerodynamic moment coefficient, about the center of the splierc is given 

as a function of angle-of-attack in Figure IV-13. 4>s expected, this configura- 
tion is statically stable at a = 0 since C = 0 and 9 / 9a < 0 at a = 0, i. e. , 

M 

if the body is disturbed from o = 0, an aerodifuaAite moment is created w hich 
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tends to drive the body back towards a = 0. However, it is not a sure thing that 
the body will come to rest again at a = 0; it might oscillate with constant ampli- 
tude (zero dynamic stability) or ^viih increasing amplitude (djoiamic instability). 
Therefore, the dynamic stability of the configuration must be examined. 

2. DYNAMIC STABILITY 

Consider angular oscillations (about center of mass) in a plane (single 
degree of freedom motion) \vith only aerodynamic torques acting, except for 
the Initial disturbance torque. The complicated dynamic interaction between 
the tether wire and the satellite is ignored. This idealized model might be use- 
ful for examining oscillatory motions in the yaw plane. 

The one-dimensional equation of rotational motion is 


la = 



a + 



where 


(lV-1) 
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or, in familiar foi’m: 
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The quantity K is found by taking the slopes of the curves in Figure IV-13 
at q; = 0. To keep K constant and equation (lV-2) linear, the majamum angle- 
of-attack must be limited to approximately 20®. The quantity 


ac 


M 


^ref ) 

u 


q=fl!=0 


is the aerodynamic damping derivative. Its value was estimated at -2.0 from 
Reference IV-6 which uses the Newtonian impact theory that is valid for con- 
tinuum flow. No attempt was made to modify the damping derivative to account 
for rarefied flow. 

The solution to equation (lV-2) is 


-f ‘ . 




a = sin 
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'Che damped natural frequency of oscillation is 

' ■ i; /f - (§)'' 

and the undamped natural frequency is 


N 
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The aerodynamic damping was so small that there was no significant difference 
in natural frequency between the damped and undamped cases. The natural 
frequency is given as a function of altitude in Figure IV-14. 'Hie differences 
in frequency for the different sized bodies are small. 
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K it is assumed that at t = 0 and a = 0, the satellite is given an impulsive 
angular speed, aj, the maximum angular deflection is given by 


“max 



(lV-6) 


The maximum angle is given as a function of the initial angular speed disturbance 
in Figure IV-15 for 107 km altitude. For a given oiq , increasing the size of the 
afterbody has an adverse effect in limiting the maximum angle; this is due to the 
increase in moment of inertia. However, for a given angular impulse disturbance. 
Figure IV-16 shows that increasing the size of the afterbody is very effective in 
limiting « 


To evaluate the effectiveness of aerodynamic damping, we find the time 
required, for the amplitude of the oscillations to decay to 1/ e of a 

From equation (lV-3) , 


D 


21 

C 


(lV-7) 


The decay time is given as a function of altitude in Figure IV-17. This time 
increases rapidly with altitude, and for altitudes greater than 110 to 115 km 
the aerod 3 mamic damping is practically ineffective. Some form of internal 
damping might be required for the higher altitudes. 

The selection of the afterbody size will depend on the disturbances 
expected, e.g, , wind profiles or jerks from the tether wire, and on the maxi- 
mum allowable amplitude of oscillation and decay time. No attempt was made 
to estimate the disturbances to be expected. 

Tlie aerodynamic heating and equilibrium surface temperature of the 
spherical satellite was given in Figures IV-5 and IV -G. The aerodynamic 
heating of the conical afterbody was obtained in a similar way from experi- 
mental data in Reference IV -2 and pi’esentcd in Figure IV-18. The data in 
Reference IV-2 are for complete cones; however, these data were used unchanged 
to estimate the convective heat transfer to the cone frustums. Performing an 
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energy balance as before yields the equilibrium temperatures in Figure IV-19. 
Separate cxirves were obtained for the different sized bodies at the lower alti- 
tudes in the transitional flow regime. 


D. Satellite Interference with Atmospheric Measurements 

If an attempt is made to measure atmospheric properties, such as 
density, by means of an instrument mounted at the surface of the satellite, 
the measurements are likely to be far different from the values in the undisturbed 
atmosphere that one would like to record. The air molecules pile up on the front 
of the satellite in a shock layer as illustrated in Figure IV-20. The density, 
temperature, and pressure are greatly Increased above their undisturbed values, 
and the air molecules are partially dissociated and ionized. One might attempt 
to overcome this satellite interference by extending a probe through the shock 
layer into the undisturbed freestream. The scope of this preliminary study is 
limited to estimating the shock layer thiclmess and, therefore, the length of 
probe required. 

A computation method for the ideal gas flow field in the stagnation region 
of a spherical body was developed by Jain and Adimurthy [IV-7] and extended 
by Jain and Kumar [IV-8] to real gas and further modified by Hendricks [IV-9] . 
The flow field along the stagnation streamline of a 1 m diameter sphere at 107 km 
altitude was computed using this program. The ideal gas version of this program 
gave reasonable agreement with experimental temperature and density profiles 
[IV-7] and the real gas version gave fairly good agreement with other theoretical 
solutions [IV-8 and IV-9], so these data can be treated with a reasonable degree 
of confidence, although a thorough comparison with experimental data is not 
available. Figure IV-21 gives the nondimensional density, temperature, and 
atomic oxygen (mass fraction) profiles along the stagnation streamline as a 
function of nondimensional distance from the surface of the sphere for the case of a 
noncatalytic wall. From this figure the shock layer extends approximately i. 1 
radii ahead of the sphere at this altitude. Therefore, a probe longer than 55 cm 
would be required to protrude into the undisturbed freestream at 107 km; longer 
probes would be required at higher altitudes. The complicated problem of the 
interference of the probe with the flow field about the satellite was not considered. 
The probe system will undoubtedly require much more study in the future. 
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Figure IV-2. Aerodynamic heating of tether mre at o: = 90° . 
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Figure IV-7. Drag coefficient of cylinders at ot = 0“ (L/D = 2) . 
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Figure IV-S. A subsatellite shape for aerodynamic stability 
( sphere-cone-frustura) . 
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Figure rV-9. Drag coefficient of sphere with cone-frustiim afterbody in orbit 

at zero angle of attack. 
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Figure IV-10. Drag force on sphere with cone-frustum afterbody at a =0°. 
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FigTire IV -11. Axial force coefficient of sphere with cone-frustum afterbody. 
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Figure IV-12. Normal force coefficient of sphere with cone-frustum afterbody. 
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Figure IV-15. Maximum amplitude of sphere-cone-frustum subsatellite oscillation 
as function of initial angular speed disturbance at a: =0°. 
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Figure IV-16. Maximum amplitude of oscillation of sphere-cone-frustum 
subsatellite as fimction of altitude. 
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Figure IV-17. Decay time for amplitude of oscillation of sphere-cone -frustum 
subsatellite as function of altitude. 
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Figure IV-20. Interference of satellite with atmospheric measurements. 




V. THERMAL STUDIES 

William P. Baker 


A. Summary 

A preliminary thermal assessment was made to determine the external 
heating rates and the resulting temperatures of a low orbit spherical satellite 
and its tether wire, and to assess the thermal protection requirements to 
effectively maintain the temperatures of the electronic packages and the antenna 
within acceptable limits. 

The resulting recommendations are to apply: (1) 3. 81 cm (1. 5 in. ) of 
high performance insulation (HPI) to the external surface of the satellite shell, 
and (2) 2. 64 cm (l in.) of high temperature reusable surface insulation 
(HRSI-LI90C) external to the satellite's antenna. 


B. Introduction 

The objectives were (1) to determine the thermal effects of aerodynamic, 
solar. Earth radiation and albedo, and internal heat upon an Atmospheric Mag- 
netospheric Plasma in Space (AMPS) tethered subsatellite and its tether wire 
and (2) to study passive thermal control techniques toward maintaining accept- 
able temperature limits. 

Other sources [V-1 through V-6J were used for information and calcu- 
tion techniques. 


C. Analysis 

The approximate aerodynamic heating rate versus altitude was detex’mined 
and plotted in Figure V-1. Figure V-2 presents a plot of aerodynamic and solar 
heat as percentiles of the total external heating versus altitude. At the design 
altitude (107 km), the satellite's exteimal heat input percentages are approxi- 
mately 97 percent aerodynamic, 2 percent solar, 0. 4 percent Kai’th I’adiation, 
and 0,0 percent albedo. These approximate percentages at 130 km are 02 per- 
cent aerodynamic, 27 percent solar, 5 percent Earth radiation, and 0 percent 
albedo. 
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These assumptions were made in computing the nominal external tem- 
perature of a spherical uncooled steel shell satellite versus orbital altitude 
(Fig. V-3). 

1. Oxidized spherical steel shell satellite 

a. Mass — 100 kg 

b. Diameter — 1.12 m 

2. Negligible convective cooling and no radiating fins — radiation from 
external surface to space. 

3. Random rotation of sphere; i. e. , no stabilizing fins for velocity 
vector control. 

4. 107 km (350 kft) nominal orbital altitude of satellite; 90 min orbital 
time, 7450 m/s orbital velocity. 

5. 50 °C maximum internal temperature limit. 

6. Transitional to free molecular flow at 90 to 110 km altitude. Free 
molecular flow at altitudes over 110 km. 

7. Equilibrium or steady state conditions. 

8. Cii’cular orbit. 


D. Results and Graphs 

An optimum thickness 3.81 cm (1.5 in. ) of high performance insulation 
could be applied external to the subsatellite shell to reduce the heat flow from 
outside to inside. 

High temperature reusable surface insulation (HRSI-L1900) can be 
applied exteiTial to the sub satellite’s antenna for effective temperature control. 
A 2. 54 cm (l in.) thickness will maintain the antenna temperature below 150°C 
at altitudes down to < 107 km. 
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Tlie greatest increase of the tether wire temperature occurs at 
altitudes below approximately 120 km. For example, if the subsatellite is 
deployed to 107 km, the tether cable material must withstand temperatures up 
to 495°C, but only 395 “C at 114 km. 

Assuming a 50°C operationr' limit internal- to the satellite (experiment/ 
instrument package), an active and/or semi-passive thermal control system 
appears to be required at the lower altitudes (130 km or less). The thermal 
control system heat load may be minimized by applying passive high per- 
formance insulation external to the experiment/instrument paclcage. Figures 
V -3 and V -4 show typical calculated temperatures for the external surface 
and the results of reducing the internal heating load by applying HPI to the 
satellite shell. The HPI retards tlie heating rate to give long thermal transients 
to the internal portions of the satellite. Figure V-5 shows the results of 
applying 3. 8 cm (l. 5 in. ) of HPI (optimum thickness) on the transient heat up 
of the experiment/instrument package. The results indicate 30 hr to achieve a 
50°C temperature ri.se at the design altitude of 107 km. Therefore, for short 
experiment times the transient nature of the passive HPI may be adequate for 
thermal control. 

Calculations were made to determine if high temperature reusable sur- 
face insulation (HRSI-L1900) can be used to limit the satellite antenna's tem- 
perature to its maximum .allowable value (150°C) and to define an optimum 
HRSI thickness. 

Figure V-(i indicates ste.ady state antenna temperatures for different 
thicknesses of the IIRSI at the design attitude (107 km) and also for 100 km. 

The antenna temperature buildup at 100 km and 107 km altitudes using 2. 54 cm 
(1 in. ) optimum HRSI thiclmess is shown on Figure V-7. Figure V-8 indicates 
the basic construction of ,a HRSI tile. 

A temperature profile of the tether wire was calculated for 100 to 220 km 
altitude (Fig. V-9), making these .assumptions: 

1. Wire's cylindrical axis at 90° to line of flight. 

2. Sol.ar angle of incidence = 90°. 

3. No heat conduction along the wire. 

4. No electrical curi’entflow in wire. 

5. Ste.ady state condition. 

G. Absorptance (o) = 0.9 and c.nutaiice (e ) =0. 8 for the wire surface. 

7. Aerodynamic and solar heat input to tlie in-ojected area (2rf ) and heat 
output by radiation from the entire cylindrical area to deep dark .space. 



E. Conclusions 

At some flight conditions, it appears that passive thermal control methods 
will suffice to successfully maintain the temperatures of the subsatellite instru- 
ments within acceptable limits. A more indepth shady nust be made to deter- 
mine applicable ranges of altitude and/or internal heat loads. 

Active cooling will be required at the design altitude, 107 km, with 
the expected orbital and instrument heat loads. 

Additional work is needed to refine the thermal analysis for uncontrolled 
satellites. The altitudes should be determined where passive thermal control 
is adequate to limit the internal temperature < 50°C. Active cooling methods 
also should be investigated. 

Passive thermal control schemes should be evaluated for velocity vector 
controlled satellites, e.g. , finned areas for thermal radiation, heat pipe applica- 
tion, and external coatings. 
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Figure V-3. Satellite external temperature versus altitude. 
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Figure V-4, High performance insulation thickness versus heat transfer 
rate for spherical tethered satellite. 
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Figure V-7. Antenna temperature versus time. 
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0.02 in. BOROSILICATE GLASS COATING, 
p - 0.15 Ib/ft^ 


1 in. SILICA INSULATION (HRSI-LI900), p- 9 Ib/ft^ 


0.01 in SILICONE ELASTOMER, p - 88 Ib/ft^ 


0.06 in. FELT STRAIN ISOLATOR PAD, P - 20 Ib/ft^ 


0.015 in. SILICONE ELASTO.MER 


TEMPERATURE CONTROLLED STRUCTURE 
(t.9., VEHICLE SKIN. ANTENNA) 


Figure V-8. HRSI (LI900) tile — basic construction- 
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Figui’e V-9. Wire temperature versus altitude. 
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VI. TETHERED SUB SATELLITE 
COA/IMUNICATIONS SYSTEM DESIGN 


J. A. Dunkin 

A. Summary 

1. ..ntenna pattern distortions due to the tether do not appear to be a 
problem if the distance between the antenna and the tether is greater than two 
wavelengths. 

2. Antennas suitable for the Atmospheric Magnetospheric Plasma in 
Space (AMPS) tethered subsatellite have been identified which can withstand 
180®C, a temperature higher than the predicted temperature of the antenna when 
covered with LI-900 insulation. 

3. The plasma cutoff frequency may preclude the use of S-band com- 
munications at altitudes below 110 km. 

4. Most equipment for S-band operation could be off-the-shelf hardware 
with a minimum of modification, 

5. Communications at frequencies above S-band will require extra sup- 
port equipment on the AMPS facility. 


B. Introduction 

The primary communication problem, which was addi'cssod during the 
AMPS tethered subsatellite feasibility study, was to determine if a communica- 
tion link could be established from the tethci’cd subsatellite to tl’.o Orbiler using 
low cost, off-the-shelf hardware. Some associated problems wliich were 
addressed during the study were the effect of the tether on the antenna pattern, 
the temperature of the subsatellite antenna, and the effect of the plasma generated 
by the subsatellite on the communication link. 

RF tracking of the subsatellite was not considered as a part of this study; 
however, use of a scanning laser radar which is curi'cntly under development is 
suggested for missions which have critical tracking requirements. 



C. Method 


The communication requirement of the tethered subsatellite are very 
similar to the communication requirements of sounding rocket programs; 
therefore, a good selection of flight proven, off-the-shelf communications hard- 
ware is available, especially in the S-band. Since the pi’oposed AMPS facility 
will have S-band communications equipment available, a subsatellite to Orbiter 
communications link margin calculation was performed for a typical S-band 
system made up of off-the-shelf sounding rocket type components. The calcula- 
tions were based upon the following key assumptions: 

1. Subsatellite antenna gain — -2 dB 

2. Receive antenna gain — +24 dB 

3. Data rate — 200 kbps 

4. Subsatellite transmitter power — 1 W 

5. Bit error rate — 10“® 

6. Receive system noise figure — 7 dB 

7. Plasma cutoff is not a problem 

8. Maximum range from Orbiter to subsatellite — 120 km 
Assuming the above, the system will have a +7.7 dB mai’gin which is adequate. 

The effect of the plasma genei’ated by the subsatellite was evaluated by 
the method outlined by Mitchell in an article entitled, "Communications-System 
Black-Out During Reentry of Large Vehicles," Proc. IEEE, Volume 55, No. 5, 
pp. G19-G26, May 19G7. It was found that the S-band communications system 
would bo blacked out at an altitude just below 105 km. Until better data are 
available on the plasma generation charactci’istics of the subsatellite, a con- 
servative estimate is that the S-band system will operate down to an altitude 
of 110 km, A graph of plasma cutoff fx-equcncy versus altitude in kilometers 
is given in Figure VI-1. 
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Selection of a specific antenna for the subsatollite cannot be made until 
the shaoe and size of the subsatellite are well defined, since the antenna and 
outer sldn of the subsatellite are integral parts of each other. However, several 
types of antennas such as planer spirals and flush mounted wrap-around antennas, 
which are well suited to this application, are available commercially. Because 
of the surface heating encountered by the subsatellite, the antenna must be 
covered with a thermal protective shield. An excellent material for such a 
shield is LI-900. LI-900 combines low dielectric constant and low loss tangent 
with good thermal insulating properties. The use of LI-900 may require a 
slight retiming of the antenna, which will not be a major problem. 

The MSFC antenna group was consulted on the effect of the tether on the 
antenna pattern. It is their best judgment based on past expex'ience, that antenna 
patteio distortions due to tlie tether wdll not be a problem if the distance between 
the antenna and the tether is greater than 2 wavelengths. 


D. Conclusion 

Communications from the tethered subsatellite to the Orbiter using off- 
the-shelf equipment is feasible at altitudes above 110 km. At altitudes below 
110 lorn, the plasma cutoff frequency will force the use of C-band for communica- 
tions, thereby, limiting the available equipment. 

Additional work should be ilouv -o refine the data on the RF losses due 
to the plasma generated by the subsatellitc, and additional work should be done 
to define a communication system operating at C-band or higher for operation 
below 110 Ivin. Also additional work needs to be done to define systems to meet 
the tracldng requirement of specific missions. 
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Figure VI-1. Plasma cutoff frequency. 


APPENDIX 

COMPUTATIONAL ROUTINE FOR ELIMINATING 
UNKNOWN LAGRANGE MULTIPLIERS 


The purpose of this appendix is to present an algoritiun that eliminates 
the unknowTi Lagrange multipliers. These multipliers enter Lagrange’ s equa- 
tions of motion by the introduction of extraneous coordinates. As a result, the 
dynamic analysis contains more coordinates than there are degrees of freedom 
of the djmamical system. This mode of operation may not seem like progress 
because it increases the number of equations to be solved, however, the ensuing 
equations are of greater simplicity and symmetry, hi addition, it enables one 
to calculate internal forces and djmamic loads of interest. 

To eliminate the Lagrange multipliers, the equations of motion are 
w'ritten in state space form. Thus, they make up a set of ordinary nonlinear 
differential equations of first order as 


Dx = Fj, 


+ 




B 


T 


X 


(A-1) 


where represents the extenial forces acting on the system and F^. the inertial 
forces. The internal reaction forces are calculated from the relationship: 




-B X 


(A-2) 


The constraint matrix B is defined by the constraint equation as 


<I> = Bx-b = 0 


(A-3) 


The Lagrange multipliers can now be eliminated by cari’ying out the following 
steps. First, the equations of motion, equation (A-l) , are prcmuitiplied by 
the inverse of the coefficient matrix D, which yields 


A-l 



X = d-‘(f^ + Fj-b'^2^ 


(A-4) 


Secondly, the equations of constraint, equation (A-3), are differentiated with 
respect to time, which yields 


Bx + Bx 



(A-5) 


Inserting equation (A-4) into equation (A-5) results in 


-I “ ^ k = 0 (A-6) 

or 

-BD“‘b^A = b - Bx - BD"‘(Fj^ + Fp (A-7) 

and finally 

A = -(BD“‘ B^)~‘ [b - Bx - BD"‘ (F^ + rpi (A-8) 


The reaction forces themselves can be rejidily obtained from equation (A-2) as 




= +B (BD~‘ B ) 


[b - Bx -BD’‘ (Fg+ Fj)l 


(A-9) 


The mathematical analysis of the tethered subsatellite docs not make use of this 
algorithm because the unknown constraint force, which is the tension in the 
tether, can be directly obtained from the tether control law and the viscoelastic 
forces of the tether. 
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